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ABSTRACT 

We compute radial color gradients produced by an outwardly propagating circular 
wave of star formation and compare our results with color gradients observed in the 
classical ring galaxy, the "Cartwheel". We invoke two independent models of star for- 
mation in the ring galaxies. The first one is the conventional density wave scenario, 
in which an intruder galaxy creates a radially propagating density wave accompanied 
by an enhanced star formation following the Schmidt's law. The second scenario is a 
pure self-propagating star formation model, in which the intruder only sets off the first 
burst of stars at the point of impact. Both models give essentially the same results. 
Systematic reddening of B — V, V — K colors towards the center, such as that observed 
in the Cartwheel, can be obtained only if the abundance of heavy elements in the star- 
forming gas is a few times below solar. The B — V and V — K color gradients observed 
in the Cartwheel can be explained as a result of mixing of stellar populations born in a 
star-forming wave propagating through a low-metallicity gaseous disk, and a pre-existing 
stellar disk of the size of the gaseous disk with color properties typical to those observed 
in nearby disk galaxies. 

Subject headings: galaxies: individual (Cartwheel) — galaxies: colors — galaxies: stellar 
content — stars: formation 
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1. INTRODUCTION 

Ring galaxies arc believed to be the result of a 
head-on galaxy-galaxy collision. Such a collision gen- 
erates an outwardly propagating wave of density en- 
hancement in the disk of the larger galaxy, referred 
to as the target. Numerical models, starting from the 
pioneering paper of Lynds & Toomre (1976), confirm 
this picture. Sinnilations reproduce some of the ba- 
sic features observed in ring galaxies, such as rings or 
crescent-shaped structures. 

The density enhancements are likely to trigger star 
formation, and a ring of star formation is expected to 
form as a response to a circular density wave. The 
star forming ring moves with time to progressively 
outer radii as the density wave propagates outwards, 
and leaves behind an evolved stellar population with 
the youngest stars located at the current position of 
the wave. The star formation history in a propagating 
wave is preserved, therefore, in the radial color distri- 
butions of stars in the disk. The inner disk is hence 
expected to be redder than its outer parts. Colors sys- 
tematically reddening towards the center were indeed 
observed in the Cartwheel galaxy (Marcum, Apple- 
ton & Higdon 1992, hereafter MAH), thus confirming 
this basic picture. More recently, similar radial color 
gradients were also found in other ring galaxies (Ap- 
pleton & Marston 1997). 

Although there seems to be a good qualitative 
agreement between the coUisional density wave model 
and the observed color gradients in ring galaxies, 
there is not yet a model of star formation that can 
quantitatively explain the observed color gradients. 
In their attempt to explain the color gradients in 
the Cartwheel, MAH had assumed a constant star 
formation scenario (Struck-Marcell & Tinsley 1978), 
and an instantaneous burst of star formation (Char- 
lot & Bruzual 1991). The agreement between these 
models and the data is far from satisfactory. Both 
the models fail to reproduce the observed sequence 
of the reddening in the Cartwheel's disk. Part of 
this mismatch is likely due to a lack of consistency 
between the data and the models they used, as dis- 
cussed below. MAH measured the color gradients in 
the Cartwheel in nine annuli with the width w 1.3 kpc 
each. A wave propagating with the velocity of 60 90 
kms~^ takes 1.4 — 2.2 x 10^ years to cross an annu- 
lus. This is clearly inconsistent with the assumption 
of the instantaneous burst model that all the stars in 
an annulus were born simultaneously. The problem 



is more serioiis if the wave propagates more slowly. 
Constant star formation models assume an extended 
duration of star formation at each ring, whereas in 
an expanding density wave the star formation stops 
once the wave sweeps across that radius. Solar metal- 
licity had been used in both the models, whereas the 
Cartwheel's mctallicity is about a tenth of solar (Fos- 
bury & Hawarden 1977). The metallicity is one of 
the basic parameters which controls stellar evolution 
(Schallcr et al. 1992) and hence color gradients are 
expected to depend on the metallicity. 

In this paper we re-address the theoretical mod- 
eling of radial color gradients in ring galaxies taking 
into account the age-spread within a sampling annu- 
lus, and the metallicity effects. We invoke two inde- 
pendent scenarios of star formation in ring galaxies. 
The first one is the conventional density wave sce- 
nario. In this scenario, the intruder's passage through 
the disk of the target galaxy creates a density wave 
which compresses the gas, thus inducing star forma- 
tion. Star formation rate at each radius in this model 
is basically controlled by the strength of the density 
wave at that radius and the gas density. We use the 
Schmidt's law (Schmidt 1959) to relate star formation 
with the gas density. In the second scenario, we in- 
voke a pure self-propagating star formation model, in 
which star formation is dissociated from the effects of 
the density wave model. An intruder has only a role 
of triggering the first event of star formation close 
to the point of impact; phenomena related to mas- 
sive stars control the star formation at aU other radii. 
As a particular case of the latter scenario we use the 
"fire in the forest" (hereafter FIF) model earlier ap- 
plied for the interpretation of the surface brightness 
and the chemical abundance gradients in ring galaxies 
(Korchagin et al. 1998; Korchagin et al. 1999). 

The kinematical data of ring galaxies seem to favor 
the density wave model for ring galaxies. However, 
the density wave model has a shortcoming in the case 
of the Cartwheel — none of the companion galaxies 
of the Cartwheel have enough mass to generate den- 
sity waves strong enough to form a giant ring such 
as that observed in the Cartwheel (Davies & Morton 
1982). New mass estimates have not improved the sit- 
uation (Higdon 1996). Under the FIF scenario, star 
formation can be sustained subsequently even if the 
density wave is weak or totally absent. FIF model 
completely neglects, however, any dynamical effects, 
which is a serious shortcoming of the theory. Both 
models, however, suit equally well to our purpose to 
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interpret optical color gradient in a radially propagat- 
ing star forming wave. 

Although we treated the two scenarios as com- 
pletely independent in this work, the two mechanisms 
may operate simultaneously. FIF phenomena can 
play a role even when density wave by itself is strong 
enough to trigger star formation. The real situation 
may be a combination of the two scenarios. 

In Section 2, we formulate the models of star for- 
mation in ring galaxies. In Section 3, we use these 
models to compute color-color diagrams for stellar 
populations born in the wave of star formation, prop- 
agating in a purely gaseous disk with different metal- 
licities. In Section 4, we discuss the comparison of the 
theoretical radial color gradients with the observed ra- 
dial color gradients in the disk of the Cartwheel. The 
results are discussed in Section 5. 

2. STAR FORMATION SCENARIOS IN RING 
GALAXIES 

2.1. Density Wave Induced Star Formation 

Classical models of ring formation through colli- 
sions between galaxies predict an outwardly propa- 
gating circular density wave in the disk of the target 
galaxy. The density wave is expected to compress the 
gas, thus increasing the gas densities above the crit- 
ical values required for star formation. We assume 
the star formation rate is proportional to the square 
of the density enhancement C{r — Vt) propagating 
outwardly in a homogeneous gaseous disk. The bal- 
ance of the surface density of stars (Af,) at a distance 
r from the center propagating with a velocity V and 
at time t can be then written as follows: 



dMs 
dt 



-D + kbC^{r-Vt), 



(1) 



with the density enhancement C given by a Gaussian 
function: 



C{r-Vt) = Aexp 



(r - Vt)^ 



(2) 



The first term D in the right-hand side of equation (1) 
gives the decrease of stellar density due to the death 
of stars, and the second term describes the growth 
of density of stars as a result of star formation. The 
coefficients b and A; < 1 determine the rate of star 
formation and its efficiency. I is the width of a star 
forming ring. 



Any dependence of the density enhancement on the 
strength of the density wave is not included in the 
above expression. Generally, the circularly expanding 
density waves are expcictcid to change their amplitude 
with time. Numerical simulations by Athanassoula et 
al. (1997) and Struck-Marcell & Lotan (1990) show 
that the relative amplitude of the density wave re- 
mains near-constant during wave propagation. The 
amplitude of the density wave approximately follows 
the surface density profile of an unperturbed disk 
which justifies our assumption of constant density 
wave amplitude for the flat surface density profile 
used in Section 3. If the pre-coUision Cartwheel was a 
normal disk galaxy with an exponential surface den- 
sity distribution, the amplitude of the density wave 
should decrease more than an order of value when the 
wave reaches the present position of the outer ring of 
the Cartwheel. 

2.2. FIF Wave of Star Formation 

The propagation of FIF star formation wave in a 
gaseous medium with an initial surface mass density, 
Mc is described by the set of equations (Korchagin et 
al. 1998): 

^ = -D+kaMc{x,t-T) J dx' f{x-x')M,{x' ,t-T), 

(3) 

^ = -aMe(x,t) |rfx'/(> 

Equations (3) and (4) give the rate of increase of sur- 
face density of the stars Ms, and the corresponding 
decrease of surface density of the star forming gas Mc 
respectively. The first term D in the right-hand side 
of equation (3) describes the decrease of stellar den- 
sity due to the death of stars, and the second term 
describes the growth of density of stars as a result of 
induced star formation. The coefficients a and fc < 1 
determine the rate of star formation and its efficiency. 
The function / models the nonlocal "influence" of the 
star complexes located at point x' on the interstellar 
medium at point x. This function has a characteris- 
tic "radius of influence" L, which, together with the 
characteristic time of the formation of stars T, con- 
trols the velocity of the star-forming wave. 

2.3. Model Parameters 

For the density wave scenario, computations were 
carried out for propagation velocities of 25, 55, 90 
and 120 kms~^. The efiiciency of star formation k 



(x-x')M,(x',i). (4) 
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was taken to bo 0.1. The amplitude A of the density 
wave was chosen to be A = 1.5 x 10'' M0/kpc^. 

The set of parameters, as well as the particular 
choice of the function / for the FIF scenario, were 
taken from Korchagin et al. (1998). Specifically, L = 
2 kpc, the delay parameter T = 1.8 x 10^ yr, the 
efficiency of star formation k = 0.1. We refer the 
reader to Korchagin et al. (1998) for more details. 

The parameters a and b determining the rate of 
star formation were chosen so as to reproduce the ob- 
served Ha surface brightness profile in the Cartwheel 
(Higdon 1995). This gives a value of a = 0.75, and 
b = 0.025 when masses, time scales and distance 
scales arc expressed in units lO^M©, 10^ yr and 1 
kpc respectively. With such values of parameters a 
and b, the theoretical rate of star formation in the 
wave is about 1.6 x 10~^ Mqyt^^ pc^^ which is com- 
parable to the estimated rate of star formation in the 
Cartwheel ~ 2 x 10"'^ Mgyr'^ pc'^ (Higdon 1995). 
Note that the colors of the star-forming waves are lu- 
minosity ratios by definition and do not depend on the 
free parameters controlling the rate of star formation. 

The prc-collisional gas density is assumed to be 
uniform for purposes of studying the general color 
properties of the star-forming waves. In Section 
4.3, wc discuss the effects of wave propagation in a 
medium with an exponentially decreasing density pro- 
file. The colors of the star forming wave propagating 
in a medium with exponential surface density profile 
do not differ significantly from those computed for the 
flat surface density profiles. 

3. COLOR GRADIENTS IN A STAR-FORMING 
WAVE 

3.1. Population Synthesis Model 

The mass Mg of stars formed per unit area is 
distributed into individual stellar masses using the 
Salpeter's IMF with a = 1.35, and the stellar mass 
interval of 0.1 Mq < < 100 Mq. Once the IMF is 
fixed, the birth and death rates of stars can be deter- 
mined at each radial grid zone from equations (l)-(4). 
The luminosity at a given band A after time t can be 
computed as : 

LA(i) =I]I]'^KT)^(m,r,i), (5) 

m r 

where N{m, r, t) is the number of stars per unit sur- 
face area at time t with mass m and age r, and 



lA{rn. t) is the luminosity of a star in a band A. Stel- 
lar luminosities /a(to,t) are obtained using stellar 
evolutionary (Schaller et al. 1992) and atmospheric 
(Kurucz 1992) models. To determine the dependence 
of the color gradients on the chemical abundances in 
the star-forming gas, the initial metallicities adopted 
ranged from Z0/2O to 2Zq. The population synthe- 
sis technique used in this work is explained in detail 
in Mayya (1995). The results of this code are com- 
pared with those from other existing codes by Chariot 
(1996). Recent updates in the code are described in 
Mayya (1997). 

To study general color properties of stellar pop- 
ulations produced by the waves of star formation, 
we modeled the propagation of the density wave in 
a purely gaseous homogeneous disk. To do this, 
we divided the disk into thirteen annuli with widths 
1.25 kpc each. The first four annuli represent the 
"nucleus", and the outer annuli correspond approxi- 
mately to the positions of annuli in Figure la of MAH. 
The approach was used both for the density wave, and 
FIF-wave models. The Equations (l)-(5) were solved 
numerically, and the average colors were computed in 
each annulus and in the "nucleus" for the different 
velocities of the wave propagation, and for the differ- 
ent metallicities of the star-forming gas. Specifically, 
we computed color gradients produced by the wave, 
propagating with the velocities 25, 55, 90 and 120 
kms~^, and for the gas metalhcities equal to Z0/2O, 
Zo/5, Z0/2.5, Zo, and 2Z0. 

3.2. Results for Density wave and FIF wave 

Figures 1-4 show the results of these simulations 
for the density wave model at four different assumed 
velocities. Each frame in Figures 1 4 presents the 
color-color diagram computed for the moment when 
the star-forming wave reaches 16 kpc, the present po- 
sition of the outer ring in the Cartwheel. 

Two conclusions can be drawn from the theoreti- 
cal color diagrams. Color properties of stellar pop- 
ulations produced by the wave of star formation do 
not depend strongly on the velocity of the wave, but 
are sensitive to the heavy element abundance of the 
star-forming gas. Stellar populations in the wake of 
the star-forming wave have spatially sequenced colors 
similar to those observed in the Cartwheel if the gas 
metallicity is sub-solar. Figures 1-4 show that color 
diagrams do not have regular reddening towards the 
center of the disk if metallicity of the star-forming 
gas is solar or higher. This result is consistent with 
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the mcasiircmcnts of the elemental abundanees in the 
Cartwheel. In this galaxy, the heavy element abun- 
dances are deficient by a factor of ten as compared to 
the Orion Nebula (Fosbury & Hawarden 1977). 

The origin of color sequentiality in a low-metallicity 
star-forming wave lies in the dependence of color- 
magnitude diagrams upon metallicity. Lowering Z 
causes individual stars to become both brighter and 
hotter. The relative contribution of supergiants to 
the total continuum considerably decreases with de- 
creasing of metallicity. E.g., the relative contribution 
of supergiants to the total continuum in R-, and V- 
bands in a 5 20 Myr starburst is about 20-80 percent 
for solar metallicity, and decreases to 3-40 percent for 
Z = Zq/10 (Mas-Hesse & Kunth 1991). These fac- 
tors, together with smoothening of colors by a non- 
coeval star formation in photometric rings, lead to the 
ordering of color-color diagrams in a low-metallicity 
gas. At solar metallicities photometric rings contain- 
ing red supergiants break the spatial ordering of col- 
ors. Uncertainties of supergiant models do not allow 
to make exact predictions, but general trend of order- 
ing of colors in a low-metallicity star-forming wave is 
robust. 

For comparison purposes, we also computed colors 

for a FIF scenario, in which all the effects of a den- 
sity wave model are completely excluded. Figure 5 
presents color diagrams for the FIF wave propagat- 
ing with the velocity 90 kms^^ computed by solving 
equations (3), (4) and (5). The IMF, stellar mass 
intervals, and the chemical abundances of the star- 
forming gas were chosen to be the same as in the 
density wave simulations shown in Figure 3. Com- 
parison of these two figures further demonstrates that 
the color properties of stellar populations born in a 
star-forming wave do not depend on the particular 
mechanism of the wave propagation. This is because 
wave velocity and the present radius of the ring basi- 
cally determine the mean age, and hence the color, of 
the stellar populations at a given radius. 

We studied the dependence of color gradients on 
the particular choice of the initial mass function. 
Namely, we calculated color gradients for the Salpeter 
IMF with stellar mass interval of 0.5 Mq < nis < 
50 Mq. Such an IMF produces more luminous stel- 
lar populations for a given stellar mass. Simulations 
show that colors are not much affected by the choice of 
the stellar mass intervals of IMF. The choice of IMF 
affects, however, luminosities of stellar populations 
born in the wave. In general, radial color gradients 



arc strongly aff(X'ted by abundance of heavy elements 
in the star-forming gas, and to less extent are affected 
by the radial velocity of wave of star formation. Col- 
ors are insensitive to a particular mechanism of the 
star-forming wave, and to the mass intervals of IMF. 

4. APPLICATION TO THE CARTWHEEL 

Theoretical colors in the wave are bluer compared 
to the colors observed in the Cartwheel galaxy (Fig. 
8 of MAH), and hence the Cartwheel's colors cannot 
be explained by a wave of star formation propagating 
in a purely gaseous disk. This implies that a detailed 
modeling is required to fully understand the color gra- 
dients in the Cartwheel. In this section, we describe 
the method that we have adopted to explain the ob- 
served colors in the Cartwheel. 

4.1. Extinction Correction 

In principle, the redder colors of the Cartwheel in 
comparison to those of the expanding wave models 
can be the result of an increasing dust content to- 
wards the nucleus. There is evidence for both cold 
and warm dust in the Cartwheel. The cold dust is 
inferred from the HST images, which show a compli- 
cated network of dust lanes in the inner ring regions. 
The warm dust is inferred in the center regions and 
in part of the outer ring by ISOCAM observations 
(Charmandaris et al. 1999). However, there has not 
been an estimation of the extinction gradient from 
the above studies. The ISOCAM emission originates 
from non-equilibrium heating of Very Small Grains 
and Polycyclic Aromatic Hydrocarbons which do not 
contribute much to the dust mass and hence to the op- 
tical extinction (Mayya & Rengarajan 1997). There- 
fore even the high central emission in ISOCAM bands 
cannot be taken as evidence for a higher dust mass at 
the center. 

We now investigate whether the observed gas den- 
sities in the Cartwheel are consistent with the as- 
sumption that the extinction plays a major role in 
creating an ordered color gradient. We hence esti- 
mated the color excess Eb-v based on the hypothesis 
that the difference between the observed color gradi- 
ent and that for the wave propagating on a purely 
gaseous disk, can be completely explained by an ex- 
tinction gradient. Details of the wave model that we 
used for the Cartwheel will be explained in Sec. 4.3. 
Using the galactic dust/gas ratio (Bohlin et al. 1978), 
we then estimated the expected total surface density 



5 



of gas 'S'jff^fj-^^ . The corresponding values of Eb—v 
for each annulus are hsted in column 2 of Table 1, 
and the estimated surface densities of gas 'E'^^f_^_fj^ 
are given in column 3. Columns 4, 5, and 6 in Table 
1 give the observed surface densities of atomic hydro- 
gen molecular hydrogen S^*, and the sum of 
atomic and molecular hydrogen in each an- 
nulus (Higdon 1996; Horellou et al. 1998). The annuh 
corresponding to the inner ring and two outer rings 
(VIII and IX) have the observed gas densities consis- 
tent with the extinction-estimated values. In the rest 
of the annuli, the extinction-estimated gas densities 
are overestimates, implying that the real extinction 
value is lower than that required to explain the ob- 
served colors. Hence the extinction gradient cannot 
be responsible for the observed color gradient. Addi- 
tional argument against large amounts of extinction is 
low metallicity of gas which decreases dust /gas ratio. 
This agrees with de Jong (1996) who foimd that the 
reddening due to dust does not play a significant role 
in the explanation of color gradients in the disks of 
spiral galaxies. The color gradients in galactic disks 
result from the combination of differences in stellar 
ages and the metallicity gradients across the disks. 

The observed gas densities, on the other hand, can 
be used to derive the visual extinction at each an- 
nuli. The galactic dust/gas ratio gives us a value of 
0.32 mag at the annuli II VII and 1.2 mag at the 
two outermost annuli. Using Balmer line ratios, Fos- 
bury & Hawarden (1977) found the extinction in two 
outer HII regions to be Ay = 2".l. Higdon (1996) de- 
rived an upper limit to visual extinction in two outer 
ffll regions CW 17 and CW 23 of Ay = 1™.6 and 
Ay = 2"*.l which are consistent with the estimate by 
Fosbury & Hawarden (1977). Unfortunately, very lit- 
tle is known about the azimuthally averaged visual 
extinction, which is the directly relevant quantity for 
determining the corrections to the radial color gra- 
dients. The values reported by Fosbury & Hawar- 
den (1977) and Higdon (1996) are greater than the 
observed visual extinction in other extragalactic HII 
regions. Kaufman et al. (1987) found an averaged 
value of < Ay >= 1^.1 ± 0.4 for 42 HII regions in 
M81. The azimuthally averaged visual extinction in 
the Cartwheel's outer ring can also be lower than 
the observed values for two most luminous HII re- 
gions. Visual extinction estimates based on the ob- 
served gas densities give also a lower value, which is 
closer to the value found by Kaufman et al (1987) 
than that estimated from Balmer lines. Note also 



that for Ay = 2"' .l, the visual extinction correction 
to the observed colors in the two outer rings of the 
Cartwheel gives unrealistically blue colors. These col- 
ors are bluer than the theoretical colors of a very 
young (1 Myr) starburst. We adopt therefore in our 
calculations, the values Ay = 1"*.2 for the two out- 
ermost photometric rings, and Ay = 0™.32 for the 
rest of the photometric rings, following the estimates 
based on the observed gas densities. 

4.2. The Wave Velocity 

Even though the model colors are not strongly af- 
fected by the wave velocity, we discuss here the veloc- 
ity adopted because of its importance for understand- 
ing the physical nature of the star-forming waves. 

There have been studies of kinematical motions in 
the Cartwheel galaxy using different techniques. Fos- 
bury & Hawarden (1977) used optical long-slit spec- 
troscopic data, Higdon (1996) used HI data and more 
recently Amram et al. (1998) used Ha Fabry-Perot 
data. A standard method of interpretation of mo- 
tions in the ring galaxies assumes that all the inter- 
nal motions are in the plane of the ring, and that the 
kinematical motions can be represented by a system- 
atic motion of the galaxy together with ring's rotation 
and its expansion. Using this method, different au- 
thors obtain rather different values for the expansion 
velocity of matter in the outer ring of the Cartwheel 
galaxy. Fosbury & Hawarden (1977) used optical 
spectroscopy of 5 HII regions and found a radial ex- 
pansion velocity of 89 ± 40 kms~^. Higdon (1996) 
finds that the outer ring is expanding with the veloc- 
ity 53±9 kms~^. Amram et al. (1998) using the best 
spatial resolution and much larger number of points 
compared to the previous studies find the velocity of 
expansion (13-17) ±2 kms~^. With a slightly dif- 
ferent approach, Amram et al. find the value of the 
matter expansion to be 30 ± 10 kms~^. This value is 
in poor agreement with the value 53 ±9 kms~^ found 
by Higdon (1996). The authors notice, however, that 
cold and ionized gas might suffer different expanding 
phases. 

It must be stressed that all these measurements de- 
termine the internal motions of cold or ionized gas in 
the wave of star formation, and do not determine the 
velocity of propagation of the wave itself. This is ex- 
actly similar to the kinematical studies of spiral den- 
sity waves, where the measured motions within the 
spiral arms themselves (i.e. the streaming motions) 
determine the amplitude of spiral density waves, and 
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do not give the pattern speed of the wave. In our 
modeUng, we therefore based the choice of the prop- 
agation velocity of the star-forming wave on the best 
fit for the surface brightness distributions in the disk 
of the Cartwheel. Korchagin et al. (1998) found that 
the observed Ha and red continuum surface bright- 
ness distributions are well reproduced for the velocity 
of the wave of about 90 kms~^. 

4.3. Color Gradients in the Cartwheel 

The color gradients that were presented in the pre- 
vious section assumed a homogeneous gaseous disk. 
The majority of galaxies have exponentially decreas- 
ing gas surface density profiles. Nevertheless, the 
above calculations adequately illustrate the basic de- 
pendence of color gradients as a function of the wave 
velocity and metallicity. To calculate the color gra- 
dients in the Cartwheel, we used an exponential gas 
surface density profile 

(Tc(r) = Aexp(-r/F) (6) 

with the parameters A = 1.5 x 10^ M0kpc~^ and 
H = 6.5 kpc which are typical for the late-type spi- 
rals (Tacconi & Young 1986). The initial metallicity 
of gas was chosen to be one-twentieth of the solar 
metallicity, close to the values observed for the outer 
ring by Fosbury & Hawarden (1977). 

FIF wave propagating in an inhomogeneous medium 
has a non-constant velocity which affects the ages of 
stellar populations behind the star-forming wave, and 
hence their color properties. The change in the wave 
velocity due to the radial variation of surface density 
is within 12% of the mean value 94 kms~^ through- 
out the disk of the Cartwheel. The resulting model 
color-color diagram is shown in Figure 6. Colors do 
not differ significantly from those computed in the 
previous section (see the top- left plot in Figure 5). 
Figure 6 also includes the extinction-corrected col- 
ors observed in the Cartwheel (filled triangles). The 
extinction corrections were based on the gas column 
density and the galactic A^/N-n ratio as explained 
in the previous subsection. Specifically, the observed 
colors for the two outer photometric rings were cor- 
rected for the extinction Ay = 1'".2, and the rest were 
corrected for the extinction value Ay = 0^.32. The 
model reproduces the sequence of points in the color- 
color diagram. However V — K colors in the wave are 
around 1 magnitude bluer than the observed colors of 
the Cartwheel galaxy. Hence, the color properties of 
the stellar populations born in a purely star-forming 



wave together with the extinction correction cannot 
explain entirely the observed colors. This is not un- 
expected, because the pre-coUisional stellar disk is 
known to have a non-negligible contribution to the 
colors of Cartwheel. 

Optical color gradients are common among "nor- 
mal" disk galaxies. Figure 7 shows B — V ysV — K 
color profiles in the disks of six nearby galaxies which 
were plotted using de Jong's (1996) observational 
data. The observed colors of spiral galaxies are even 
redder than the observed colors in the Cartwheel. 
This indicates that the Cartwheel's colors can indeed 
be obtained by mixing stellar populations born in the 
wave with the older stellar populations existing in the 
Cartwheel's disk before the collision. 

Color mixing cannot be considered separate from 
the surface brightness constraints. The spatial distri- 
bution of the young stellar population can be found 
from the best fit to surface brightness distribu- 
tion in the Cartwheel. Assuming exponential gas sur- 
face density distribution we find that the best fit to 
Ha surface brightness distribution in the Cartwheel 
(Figure 8) can be obtained for the FIF-wave propa- 
gating in the gaseous disk with the parameters A = 
5 X lO'^M© kpc~^, and H = llkpc. For the old stellar 
population in the Cartwheel we assumed that it has 
an exponential distribution in the F-band. Specif- 
ically, we assumed that V-band surface brightness 
/i(r) at a distance r from the galactic center is given 
by the relation /Lt(r) = iiq + 1.086r/i?0j where /Lto 
is the surface brightness in the center of the disk 
in magarcsec"^, and Rq determines the scale length 
of the exponential surface brightness distribution. 
Choosing Rq = 5 kpc which is the typical value for 
disk galaxies (de Jong 1996) we find that the best fit 
to the red continuum profile in the Cartwheel can be 
obtained for /xq = 23.5 magarcsec"^. The red con- 
tinuum brightness distribution is however in a poor 
agreement with observations. Maximum of the red 
continuum peak associated with the outer ring of the 
Cartwheel is less prominent as compared to obser- 
vations. Higher gas metallicities increase a relative 
contribution of supergiants to the total continuum. 
Figure 9 shows Ha and red continuum surface bright- 
ness profiles for the same parameters as in Figure 8 
but for the metallicity of gas Zq/5, which gives more 
satisfactory fit to observational data. 

The surface brightness profiles were then used to 
calculate the luminosity of pre-coUisional stellar disk 
of each radial zone. The total luminosity of each zone 
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is obtained by summing the luminosities of pre- and 
post- collisional stars, which were then used to calcu- 
late the color profiles. The pre-collisional stellar disk 
seems to lack any bulge, and hence, the Cartwheel was 
most probably a late- type spiral (Higdon 1996). We 
assume that the pre-existing old disk of the Cartwheel 
has color properties typical to those in the late-type 
disk galaxies. As an example we used de Jong's (1996) 
color profiles for the Sbc galaxy UGC 01305. The 
resulting color profiles are plotted in Figures 10-11 
for Z0/2O and Zq/S metallicities respectively. The 
filled triangles in the Figures indicate the extinc- 
tion corrected observed colors of the nine rings in 
the Cartwheel's disk, and its nucleus. The colors of 
galaxy UGC 01305 are marked by open circles. The 
resulting colors of the mixed stellar populations are 
shown by the open squares. The outer three radial 
zones, 9, 8 and 7 are identified both for the observed 
and model points. It can be seen that the observed 
colors are best reproduced in the model with expo- 
nentially decreasing gas surface density distribution 
and with metalicity of gas Z0/2O (Figure 10). This 
model, however, does not reproduce well the red con- 
tinuum surface brightness profile. The model with the 
metallicity of the star forming gas Z0 /5 satisfactory 
reproduces the range of the observed colors and their 
spatial sequence. This is not the case if the metal- 
licity of the stars formed in the wave was solar, as 
illustrated in Figure 12. All the notations of the pre- 
vious figure are retained, and hence the open squares 
represent the mixed model colors. This model neither 
explains the range of observed colors nor the spatial 
sequence of the colors. 

5. SUMMARY AND DISCUSSION 

This paper has focused on the optical properties 
of disk galaxies with large-scale star forming rings. 
We considered an improved model of a star forming 
wave, and used the population synthesis approach to 
calculate the disk color gradients produced by an out- 
wardly propagating wave of star formation. The re- 
sults of our calculations can be summarized as follows: 

1. Sub-solar metallicities (at least by a factor of 
2.5) are necessary for generating color profiles which 
redden systematically towards the center. This re- 
sult is consistent with the observed metallicity in the 
Cartwheel, where such color gradients are observed. 

2. The theoretical colors of the stellar populations 
born in the wave propagating in a purely gaseous 



disk are much bluer than the observc;d colors of the 
Cartwheel's disk. A reddening gradient caused by 
interstellar dust cannot explain the discrepancy. In- 
stead, an explanation of the Cartwheel's colors re- 
quires the presence of old stellar disk extending out 
to the present position of the ring. The existence of 
a pre-collisional nucleus in the Cartwheel has been 
known for more than two decades, but this is the first 
unambiguous evidence for a large stellar disk in the 
Cartwheel. 

Simple models used in this paper leave aside many 
important aspects of the dynamics of ring galax- 
ies. They allow, however, to make a conchision on 
the nature of pre-collisional Cartwheel galaxy. Hig- 
don (1996) assumed that a pre-collisional Cartwheel 
would have appearc;(l as a small spiral embedded in 
an extensive low surface density HI disk. Our results 
show that the old stellar populations must exist in 
the entire disk of the Cartwheel in order to under- 
stand the observed color gradients. 

These conclusions raise the question of the sur- 
vival of the pre-collisional disk structure during an 
encounter. Athanassoula, Puerari & Bosma (1997) 
numerically studied collisions of intruders with barred 
disk galaxies. They found that an intruder with mass 
equal to ten percent of the mass of a target barred 
galaxy leaves the bar un-destroyed. To our knowl- 
edge, there is no study of the collisional disruption 
of disk galaxies that have well-developed spiral struc- 
ture. The stability properties of gravitating disks can 
be drastically changed by a few percent admixture of 
cold gas. If the pre-collisional disk of Cartwheel was 
gas rich, and most of the gas was burnt in the star- 
forming wave or stripped from the galactic plane, any 
pre-collisional spiral structure would not have sur- 
vived. All these questions remain for future inves- 
tigations. 

The authors thank an anonymous referee for his 
useful comments, W. Wall for reading the manuscript 
which resulted in improved presentation, and Roelof 
de Jong for providing his observational data in elec- 
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his hospitality, and the National Astronomical Obser- 
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search grant 211290-5-25869E. 



8 



REFERENCES 

Amrani, P., Mcndcs dc Oliveira, C, Boulesteix, J., & 
Balkowski, C. 1998, A&A, 330, 881 

Appleton, P.N., & Marston, A.P. 1997, AJ, 113, 201 

Athanassoula, E., Puerari, I., & Bosma, A. 1997, MN- 

RAS, 286, 284 

Bohlin, R.C., Savage, B.D., & Drake, J.F. 1978, ApJ, 
224, 225 

Chariot, S. 1996, in ASP Conf. Ser. 98, From Stars 
to Galaxies, Ed. C. Leitherer, U. Fritze-v. Al- 
venslebcn, & J. Huchra (San Francisco: ASP), 275 

Chariot, S., & Bruzual, A. 1991, ApJ, 367, 126 

Charmandaris, V., Laurent, O, Mirabel, I.F., Gallais, 
P., Sauvage, M., Vigroux, L., Cesarsky, C, & Ap- 
pleton, P.N. 1999, A&A, 341, 69 

Davies, R.L., & Morton, D.C. 1982, MNRAS, 201, 
69P 

de Jong, R.S., 1996, A&A, 313, 377 

Fosbury, R.A.E., & Hawarden, T.G. 1977, MNRAS, 
178, 473 

Higdon, J.L. 1995, ApJ, 455, 524 

Higdon, J.L. 1996, ApJ, 467, 241 

Horellou, C, Charmandaris, V., Combes, F., Apple- 
ton, P.N., CasoU, F., & Mirabel, I.F. 1998, A&A., 
340, L51 

Kaufman, M., Bash, F.N., Kennicutt, R.C., & Hodge, 
P.W. 1987, ApJ, 319, 61 

Korchagin, V., Mayya, Y.D., & Vorobyov, E.I., Kem- 
bhavi, A.K. 1998, ApJ, 495, 757 

Korchagin, V., Vorobyov, E.I., & Mayya, Y.D. 1999, 
ApJ, 522, 767 

Kurucz, R.L. 1992, lAU Symp. 149, Stellar Popula- 
tions of Galaxies, Ed B. Barbuy & A. Renzini (New 
York: Kluwer), 225 

Lynds, R., & Toomre, A. 1976, ApJ, 209, 382 

Marcum, P.M., Appleton, P.N., & Higdon, J.L. 1992 
ApJ, 399, 57 (MAH) 

Mas-Hesse, J.M., & Kunth, D. 1991, AAS, 88, 399 



Mayya, Y.D. 1995, AJ, 109, 2503 

Mayya, Y.D. 1997, ApJL, 482, L149 

Mayya, Y.D., & Rengarajan, 1997, AJ, 114, 946 

Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 
1992, A&AS, 96, 269 

Schmidt, M. 1959, ApJ, 129, 243 

Struck-Marcell, C, & Tinsley, B.M. 1978, ApJ, 221, 
562 

Struck-Marcell, C, & Lotan, P. 1990, ApJ, 358, 99 
Tacconi, L.J. & Young, J.S. 1986, ApJ, 308, 600 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v4.0. 



9 



V-K 



V-K 



0.4 
0.3 
0.2 

> 0.1 
0.0 

-0.1 
-0.2 

0.3 
0.2 

> ° '' 
CO 

0.0 
-0.1 
-0.2 



: v=25km/s; z=z0/20 :: v=25km/s; z=z0/5 



I I I I I I I I 1 I 



i: v=25km/s; z=zO 



i_v=25km/s; z=z0/2.5 

: 

2-1012 -1012 

0.4 H 




0.3 
0.2 'r 
^ 0.1 - 

CQ 

0.0 r 
-0.1 

-0.2 'r 



v=25km/s; z=2*zO 



-2-10123 
V-K 



Fig. 1. — The B — V, V — K color plots compute^ 
for nine annular zones of a model ring galaxy. Th 
star-forming density wave propagates with a veloc 
ity 25 kms"-^ in a purely homogeneous gaseous disl 
The plot illustrates the radial dependence of color 
within the disk for five different metallicities of the 
star-forming gas Z0/2O, Z0/5, Z0/2.5, Z©, and 2^0. 
Colors have regular reddening towards the disk center 
if the initial metallicity of gas is low. 
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Fig. 2. — Same as in Fig. 1, but for a wave velocity 
of v=55 kms~^. 
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Fig. 3. — Same as in Fig. 1, but for a wave velocity 
of 90 kms-^ 



Fig. 4. — Same as in Fig. 1, but for a fast star- forming 
wave (v=120 kms~^). 
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Fig. 5. — The B — V, V — K color plots computed for 
the FIF wave propagating with a velocity 90 kms~^. 
The color diagrams for the FIF wave are similar to the 
color diagrams of the density wave propagating with 
the same velocity (Fig 3). The top-left frame shows 
also color diagram computed for the FIF wave propa- 
gating in the disk with an exponential surface density 
distribution (open squares) used for the modeling of 
the color properties of the Cartwheel galaxy. 



Fig. 6. — The model ring galaxy colors compared with 
the observed colors of the Cartwheel. The large filled 
triangles mark the B — V , V—K colors observed in the 
Cartwheel galaxy (MAH). The colors of the two outer 
photometric rings, VIII and IX were corrected for ex- 
tinction using Ay = 1'".2. The colors of the inner 
photometric rings and the Cartwheel's nucleus were 
corrected using Ay = 0'".32. The filled circles show 
the B — V , V — K colors produced by the star-forming 
wave propagating in a purely gaseous (Z0/2O) disk 
with the velocity 90 kms~^. The model colors are 
bluer than the observed colors in the Cartwheel. The 
values of extinction correction are shown by the ar- 
rows. 



12 



V-K V-K 




^ ^ F ^ 0.5 

2.0 2.5 3.0 3.5 2.5 3.0 3.5 4.0 



Fig. 8. — The red continuum and Ha surface bright- 
ness distributions produced by the wave of star for- 
mation propagating in a gaseous disk with expo- 
Fig. 7. — The color-color {B — V, V — K) plots of six nential surface density profile with the parameters 
disk galaxies as measured by de Jong (1996). A = 5 x lO^'MQkpc^^, H = llkpc and gas metal- 

licity Z0/2O. An underlying old stellar disk with 
Mo(^) = 23.5 magarcsec"^ and i?o = 5 kpc is added 
to the model ring galaxy. 
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Fig. 9. — The red continuum and Ha surface bright- 
nesses in the wave with the same parameters em- 
ployed in Figure 8 but with metallicity Z0/5. 



Fig. 10. — The model colors computed for the wave 
of star formation propagating in a gaseous disk with 
the parameters employed in Figure 8. An underly- 
ing disk added to the model ring galaxy. The large 
filled triangles mark the extinction corrected colors 
of the Cartwheel galaxy. The open circles indicate 
the observed colors for galaxy UGC 01305. The open 
squares show the colors obtained as a result of mixing 
of the colors produced by the star-forming wave prop- 
agating in a purely gaseous (Zq/20) exponential disk, 
with the observed colors in UGC 01305. The mixed 
colors are close to the observed ones in the Cartwheel. 



14 




IX VIII 



1 2 
V - K (mag) 



1.2 
1.0 
0.8 

3 0.6 

E. 

> 0.4 
0.2 
0.0 
-0.2 
-0.4 



nucleus (UGC01305) 





IX 



VIII 



1 2 
V - K(mag) 



nucleus 
» 



Fig. 11. — Same as in Fig. 10, but for the initial 
metallicity of the gaseous disk Zq/S. 



Fig. 12. — The model colors computed for the model 
with the parameters employed in Figures 10 and 11, 
but for the initial metallicity of the gaseous disk 
Z— Z0. The observed colors are in disagreement with 
the Cartwheel's observations, suggesting a low metal- 
licity for the star-forming gas in the Cartwheel. 
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Table 1: Gas surface densities estimated from the 
maximum allowed color excesses compared with the 
observed gas surface densities. 



